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Recently, random matrix theory predictions for the distribution of low-lying Dirac operator eigen- 
values have been extended to include lattice effects for both staggered and Wilson fermions. We 
computed low-lying eigenvalues for the Hermitian Wilson-Dirac operator and for improved stag- 
gered fermions on several quenched ensembles with size « 1.5 fm. Comparisons to the expecta- 
tions from RMT with lattice effects included are made. Wilson RMT describes our Wilson data 
nicely. For improved staggered fermions we find strong indications that taste breaking effects on 
the low-lying spectrum disappear in the continuum limit, as expected from staggered RMT. 
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1. Chiral perturbation theory including lattice artifacts 

In recent years, chiral perturbation theory (XPT) has been extended to include lattice artifacts 
for both staggered [JTJ] and Wilson fermions [ph. In this talk, we will be concerned only with the 
£-regime of ZPT where the zero momentum modes dominate. Thus, we consider only the zero 
momentum part of the chiral Lagrangian 

Jgf = --mLTr(U + U f ) +a 2 Y . (1.1) 

Here ~f describes the lattice artifacts. For staggered fermions these are dominated by taste breaking 
terms [jl|] 

y = - l -C A Tv{^ 5 U^ + h.c)+... , (1.2) 

where we displayed only the term that dominates the pseudoscalar mass splittings explicitly. Here 
%H = are taste matrices. This is also the term that dominates the lattice effects in the low-lying 
Dirac spectrum in the regime of weak taste breaking. 
For Wilson fermions, the lattice artifact terms are [§] 

r = W/ 8 Tr([/ 2 + [/ t2 )+W 6 [Tr([/ + [/ t )] 2 + W 7 [Tr([/-[/ t )] 2 . (1.3) 

The two-trace terms, with coefficients We and Wj, are suppressed at large N c . We will neglect them 
here and only keep the one-trace correction term, proportional to Wg. 



2. Random matrix theory including lattice artifacts 

The £ -regime of XVI at leading order, can equivalently be described by a chiral random matrix 
theory (RMT). For continuum QCD, the Dirac operator is represented in RMT as 



iW 
iW f 



(2.1) 



with W a random (N + v) x N complex matrix, when working in a sector with index (topological 
charge) V. 

2.1 Staggered RMT 

For staggered fermions, the Dirac operator in staggered RMT (SRMT) is represented as ^ 

%ta g = %®U + a<?, '^ = {^0 B^j® ^ 5 ■ (2 - 2) 
Here I4 is the identity matrix in taste space and 3? denotes the taste beaking terms. The dominant 



taste-breaking term, corresponding to the term explicitly shown in Eq. (1.2), denoted by j5q is 
shown in the second part. There, and are random Hermitian matrices of size (A^ + v) x (A^ + 
v) and N xN, respectively. The width of the Gaussion distributions of and B^ is proportinal to 



C4. The dimensionless combination a 2 QV controls the strength of the taste breaking in SRMT. 
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Without the taste-breaking terms (i.e., in the continuum limit) the eigenvalues of the Dirac 
operator come in degenerate quartets. For weak taste breaking, a C4V <C 1, the quartets of eigen- 
values are split, at leading order in a perturbation by a^ 4 , into pairs of eigenvalues with splitting 
AX quar t [f||. The pairs of eigenvalues, in turn, are split, for weak taste breaking, at second order in 
a perturbation by a<5c 4 , with subdominant splitting AA pai >. So SRMT predicts 

^Icay^, ^2±o, a 2 C,V. (2.3) 
A A 

2.2 Wilson RMT 

For Wilson fermions, the Dirac operator in Wilson RMT (WRMT) is represented as [H] 



b + 5 [ n „ , (2.4) 




with A and B random Hermitian matrices of size (N + v) x (N + v) and N x N, respectively, that 
represent the chiral symmetry breaking Wilson term of the lattice Wilson-Dirac operator. 

Akemann et al. , Ref. [Q], have worked out the eigenvalue distribution of the Hermitian Wilson- 
Dirac operator, H\y = 75 (Dw + mo), or its RMT equivalent, = 75 + m ) with 

m = mLV = 2mN and a 2 = a 2 W 8 V = \a 2 N (2.5) 



held fixed, using Wilson ZPT The results were reproduced directly from WRMT in [g]. In Eq. ( pL5| ) 
m is a suitably subtracted version of the bare mass mo. In the analytical predictions, the eigenvalues 
are rescaled with LV in the lattice QCD case, and with 2N for WRMT. 



3. Results for staggered fermions compared to SRMT 

We have computed low lying staggered eigenvalues for various staggered actions with different 
smearings using ensembles of pure gauge configurations generated with the Iwasaki gauge action. 
As observed previously [Qj, the quartet structure, and the appearance of clearly distinct would-be 
zeromodes, becomes more visible with increased smearing and with smaller lattice spacing. Here 
we show and discuss some results with our best staggered action, a HISQ-like action but with first 
a HYP(ii) smearing step |7|] instead of a fat7 smearing step and subsequent unitarization, followed 
by an asq step with the appropriate Lepage term, denoted by "HYPasq2L". For more details and 
further results we refer to Ref. [§]. Here we illustrate a few of our findings. 

We compare distributions of the lowest four eigenvalues from ~ 1200 Q = configurations at 
a = 0.093 fm and L = 1.5 fm with Monte Carlo generated SRMT v = eigenvalue distributions 
for the taste violating parameter a 2 C<\V = 0.005 in Fig. [Tj 1 The LQCD eigenvalues are rescaled by 
LV and the SRMT ones by 2Af for the comparison. The distributions agree quite nicely. 

Fig. ^ shows the distribution of the lowest four eigenvalues, including the two positive would- 
be zeromodes, for ~ 2000 \Q\ = 1 configurations of the same gauge ensemble and v = 1 SRMT 
MC generated eigenvalue distributions. Again, the distributions agree quite nicely, although the 

1 We thank James Osborn for providing the eigenvalues from SRMT generated by MC. 
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Figure 1: Distributions of the lowest four (rescaled) eigenvalues for Q — 0. The left panel shows the 
staggered Dirac eigenvalues from an ensemble of a = 0.093 fm and L = 1.5 fm configurations, the right 
panel the eigenvalues from a MC generated SRMT ensemble with taste-breaking parameter a 2 C^V = 0.005. 



Ev_HYPasq2L, /3 Iw =2.635, 16 4 , Q=l 



Ev_$RMT, a 2 C 4 V = 0.005, v=\ 



1 



2 j;i 



x ev 1 
□ ev 2 
+ ev 3 

h ev 4 




0B- 



x ev 1 
□ ev 2 
+ ev 3 

k ev 4 



Figure 2: Same as Fig. [j] but for \Q\ = 1 (V = 1) configurations. Here, the lowest two eigenvalues are the 
would-be zeromodes. 



distributions of the would-be zeromodes for the staggered Dirac operator have somewhat longer 
tails than their SRMT counterparts. In the SRMT distributions, the peak height of the "would-be 
zeromodes" turns out to be the feature most sensitive to the value of the taste-breaking parameter 
a 2 C<\V . It is from this peak height that we estimated that a 2 C^V = 0.005 provides the best match to 
the LQCD data. 

According to the SRMT prediction, Eq. (2.3), the dominant splitting AX quart should vanish, 
at fixed volume, like a"a, where n depends on the degree of improvement, while the subdominat 
splitting AXp a i r should vanish like a 2 . We compare these predictions, as well as the vanishing of the 
would-be zeromodes, for fixed lattice size L = 1.5 fm, in Fig. |3[ The would-be zeromodes appear 
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HYPasq2L splittings for L=1.5 fm 
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Figure 3: The vanishing of the splitting (left panels) and would-be zeromodes (right panels) - shown is the 
average over the would-be zeromodes with slightly positive eigenvalues - as a function of a 2 (top) or aja 
(bottom). The inclined straight lines show the behavior to be tested in the plots. 



to vanish a little faster than a 2 . The behavior of the splittings is less clear, with a vanishing as a 2 
favored for the smaller splitting, AX pa i r . But both splittings certainly vanish in the continuum limit, 
as expected. 



4. Results for Wilson fermions compared to WRMT 2 

We computed the low-lying eigenvalues of the Hermitian Wilson-Dirac operator H w on the 
<2 = and \Q\ = 1 configurations of an ensemble with lattice spacing a = 0.075 fm and size L = 1.5 
fm, as before generated with the Iwasaki gauge action, for two bare quark masses. To suppress 
dislocations we smoothened the gauge field with one HYP smearing step before constructing the 
Wilson-Dirac operator. 

2 The results in this section were obtained in collaboration with Poul Damgaard and Kim Splittorff. 
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Figure 4: The eigenvalues density of H w with bare mass amo — —0.184 on an ensemble of a = 0.075 fm and 
L = 1.5 fm configurations with Q = Q (left) and \Q\ — 1 (right). The Q = distribution was used to obtain 
the WRMT parameters. The same parameters are used for the \Q\ = 1 WRMT distribution on the right. 



RMT predictions are made for sectors with a given index v. Such an index can be defined 
from the real eigenvalues of the Wilson-Dirac operator and the chiralities of the corresponding 
eigenmodes Ref. PJ]. Equivalently, the index can be obtained from the spectral flow of the Hermi- 
tian Wilson-Dirac operator [ph. In both cases the definition depends on a cut: the maximum real 
eigenvalue kept or the mass at which the flow is terminated. For our setup we found little depen- 
dence on this cut and good agreement with the cheaper definition of topological charge with six 



HYP smearings and use of an improved F*F operator [10]. 



To compare the sepctrum of the Hermitian Wilson-Dirac operator with predictions from WRMT, 
we used the distribution in the Q = sector with the bare mass corresponding to a lighter quark 
mass to find the WRMT parameters m and a, and the eigenvalue rescaling factor LV that best 
describe our data, see Fig. ^| (left). Using the same parameter values, WRMT then predicts the 
\Q\ = 1 distribution that can be compared with the measured one, see Fig. § (right). We find a nice 
agreement with the measured (histogrammed) distribution. 
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Figure 5: Same as Fig. [| but for bare mass amo = —0. 178. The same WRMT parameters, up to changing m 
according to the bare mass difference, 8m — 8moT.V, were used. 



Changing only the bare quark mass in H\y should leave the WRMT parameter a, and the 
rescaling factor LV unchanged, while the WRMT parameter m should be changed by 8m = Sm^LY, 
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Eq. (2.5), where 8mo is the difference in bare quark mass. We, therefore, have parameter free 
predicitions for the eigenvalue distributions with the second quark mass. As shown in Fig. ||, these 
predictions work well. 




X X 
Figure 6: The dependence of the distribuion of the eigenvalue density of Jtffc- on m (left) and on a (right). 

We illustrate the sensitivity of the analytic WRMT distributions to the WRMT parameters m 
in Fig.|(left)andto a in Fig. |6] (right). For the ensemble considered in Figs. D and || we can obtain 
m to an accuracy of about 0.5 and a to about 0.05. Going to smaller lattice spacing will increase 
the accuracy. For further details and more tests, see [O]. A. Deuzeman presented results from a 
similar study at this conference with compatible results [|l2|]. 
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